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ABSTRACT 

The  optical  absorption  profile  of  an  atom-surface  vibrational  bond  is 
studied.  From  Markovian  relaxation  theory  it  follows  that  the  line  shape  is 
a  Lorentzian  around  the  adsorbate  resonance  frequency  u^.  Dispersion  relations 
for  crystals  have  a  fairly  small  cutoff  frequency  w  ,  which  prohibits  the  use 
of  such  a  simple  theory.  We  calculated  the  spectral  profile  with  finite  memory- 
time  reservoir  theory,  and  we  found  that  the  modified  Lorentzian  vanishes  above 
Up.  Also  a  new  spectral  line  at  u^+Up  is  predicted,  which  disappears  in  the 
Markovian  limit.  The  physical  origin  of  the  new  line  is  explained. 
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Adsorbed  atoms  on  the  surface  of  a  crystal  can  absorb  photons  from  an 
incident  infrared  laser  beam.  Coupling  between  the  van  der  Waals  bond  and  the 
radiation  is  brought  about  by  the  motion- induced  dipole  moment  fi  (the  atom  itself 
is  assumed  to  be  neutral  for  IR  light),  and  the  atom-crystal  interaction  is 
governed  by  the  single-phonon  coupling  Hamiltonian1- 


Hi  “  -(u-ez)dV/dz  , 


in  terms  of  the  derivative  of  the  binding  potential  well  V(z) .  Here  z  denotes 

the  normal  direction  to  the  surface  (we  neglect  lateral  motion)  and  u  is  the 

2  ~ 

phonon- field  amplitude  operator,  evaluated  at  the  position  of  the  surface  atom 
which  is  closest  to  the  adsorbate.  Relaxation  of  the  atomic  bond  due  to  phonon 
exchange  with  the  crystal  is  most  conveniently  described  with  reservoir  theory, 
where  the  crystal  is  regarded  as  a  thermal  bath.  Then  an  equation  of  motion  for 
the  reduced  adsorbate  density  operator  can  be  derived,  in  which  the  properties  of 
the  heat  bath  only  enter  parametrically,  and  can  be  expressed  in  terms  of  the 
Fourier-Laplace  transform  of  the  reservoir  correlation  function 


2Mm_ 

JQ 


I 


<(u(r)«ez)(u*ez)>  . 


Equation  (2)  defines  the  dimensionless  function  g(z)  in  terms  of  the  mass  M  of  a 
bulk  atom  and  the  Debye  frequency  Then  it  ~an  be  shown  that  g(z)  depends 
only  on  the  dimensionless  temperature  7  -  )tf«g/kgT,  with  kg  Boltzmann's  constant. 
An  explicit  expression  for  g(r;  car  be  found  ir.  Tcf .  '■  ■ 


In  the  most  simple  theory  of  relaxation  one  adopts  the  Markov  and  secular 
approximations,  which  yields  for  the  low- intensity  absorption  line  shape  as  a 


function  of  the  laser  frequency  w 


I(w)  -  wQ(n^-n2)Re 


_ i _ 

oj-w_+ia(g(w  )+g*(  -w  ) ) 
o  o 


(3) 


This  profile  is  a  Lorentzian  with  its  peak  value  at  u  -  +  aIm(g(wQ)  + 

g *(-£>  )),  and  a  half -width  at  half -maximum  equal  to  aRe(g(o>  )  +  ))•  The 

0  o  O  v 

notation  is  u>  -  u0  “  unperturbed  adsorbate  resonance,  and  n^  —  population 

of  the  i-th  level  (n^  >  nj).  The  temperature  dependence  is  incorporated  in  g(z), 
and  the  phonon- coupling  strength  is  measured  by 


(4) 


Most  crucial  for  the  derivation  of  the  line  profile  (3)  is  the  conjecture 
that  the  phonon -amplitude  correlation  function  <(u(r) *ez) (u*e  )>  decays  to  zero 
very  fast  for  r  >  0,  which  would  justify  the  Markov  approximation.  It  is  an 
essential  feature  of  a  crystal,  however,  that  its  dispersion  relation  has  a 
finite  cutoff  frequency  which  sets  the  lower  limit  of  1/w^  on  the  decay  time. 
In  fact,  the  structure  of  the  amplitude  operator  u  is  such  that  the  correlation 
function  decays  only  as  1/r  for  r  -+  «,  rather  than  exponentially,  which  makes  a 
Markov  approximation  (zero  memory  time)  at  least  doubtful.  In  order  to 
investigate  this  problem  quantitively ,  we  have  developed  a  finite  memory-time 
relaxation  theory, ^  with  no  restriction  on  the  time  scales.  Subsequently,  we 
have  applied  this  theory  to  the  evaluation  of  the  absorption  line  shape,  with  the 
formal  result^ 


•J*  K~  V,  V.  ^  V 


I(w)  "  "#*2i  ReTrLx  W-L  +ir (w)  <ViT(w))p 


This  expression  holds  for  any  configuration  of  adsorbate  levels  and  includes 

multiphonon  transitions  as  well.  In  Eq.  (5),  p  equals  the  adsorbate  density 

operator  in  thermal  equilibrium,  and  the  Liouvillians  and  L^.  are  defined  by 

L^o  -  pa ,  Ly£7  -  [ft, a]  in  terms  of  their  action  on  an  arbitrary  adsorbate  density 

operator  a.  The  relaxation  operator  r(o>)  is  frequency  dependent,  which  reflects 

the  memory  in  the  time  regression  of  two-time  quantum  correlation  functions.  In 

the  time  domain  this  would  be  a  memory  kernel,  and  in  the  Markov  approximation 

T(u)  assumes  an  w- independent  value.  Furthermore,  a  Liouville  operator  T(w) 

appears,  which  accounts  for  the  fact  that  the  density  operator  of  the  entire 

system  does  not  factorize  as  Pa(jSOrbate  x  ^crystal  t^lerma^  equilibrium  if  the 

memory  time  is  finite.  We  call  this  term  the  initial  correlation  contribution  to 

the  line  shape,  in  contrast  to  the  term  1^.,  which  is  referred  to  as  the 

regression  part.  In  a  zero  memory- time  approximation  the  operator  T(w)  vanishes 

identically,  since  then  a  factorization  of  the  density  operator  can  be 

justified.  Both  T(w)  and  T(w)  are  complicated  operators.  Elsewhere^  we  have 

studied  the  modifications  of  the  Lorentzian  profile  (3)  due  to  the  memory 

effects,  and  we  have  derived  an  explicit  expression  for  I(o>),  pertaining  to  the 

situation  where  w  is  smaller  than 

o  D 

In  this  Letter  we  report  a  peculiar  feature  of  the  absorption  profile  (5) 
which  arises  as  a  consequence  of  the  memory  in  the  interaction.  Since  Eq.  (5) 
applies  to  any  configuration,  we  can  consider  the  case  where  the  vibrational 
resonance  is  larger  than  w^.  Then  we  scan  the  laser  over  this  resonance, 

implying  that  also  w  >  For  frequencies  larger  than  a  Markov  approximation 

can  never  be  justified,  and  it  is  imperative  to  include  the  memory  in  the  time 
evolution.  At  this  stage  we  neglect  the  small  contribution  of  multiphonon 


processes,  in  order  to  emphasize  the  significance  of  the  predicted  effects.  Then 
it  can  be  shown^  that  the  regression  part  (-  L^.)  ,  which  would  lead  to  the 
Lorentzian  in  the  Markov  approximation,  is  identically  zero  for  w  >  Hence 
any  absorption  is  due  to  the  initial-correlation  term  T(u) ,  and  therefore  a  mere 
consequence  of  the  finiteness  of  the  memory  time.  We  obtain  for  the  absorption 
spectrum 


(w)  -  a(  e 


for  cj  >  Here  we  used  that  the  population  of  the  upper  state  must  be  zero  for 

>  <jq  (so  -1),  and  the  general  property  Re  g(z)  -  0  for  jz|  >  1  of  the 

reservoir  correlation  function.  A  characteristic  example  of  I(w)  is  plotted  in 
Fig.  1. 

We  observe  that  I(w)  has  a  spectral  line  for  w  >  Recalling  that  our 

result  (6)  only  includes  single-phonon  transitions,^  this  might  appear  to  be  a 

contradiction.  We  see,  however,  that  the  line  is  not  brought  about  by  the 

transition  dipole  moment  as  usua^->  t>ut  by  permanent  dipole  moments  ^2 

and  of  the  two  levels.  From  the  diagrams  of  Fig.  2  we  then  anticipate  that 

it  could  be  possible  to  absorb  a  photon  with  w  >  and  in  such  a  way  chat 

energy  conservation  can  be  met  with  a  single-phonon  process.  The  mismatch 

between  w  and  is  simply  the  phonon  frequency.  Although  level  | 2>  cannot  be 

populated  because  of  w  >  «_,  it  can  nevertheless  be  reached  via  a  photon-phonon 

o  D 

process.  This  mechanism  gives  rise  to  the  line  at  <j  > 

Figure  2  suggests  that  the  processes  responsible  for  this  line  (b  and  c)  can 
only  occur  for  |w-w  |  <  a>^,  since  otherwise  the  phonon  frequency  would  become 


larger  than  This  can  be  shown  explicitly  from  Eq.  (6).  If  we  work  out  the 

real  part  and  use  Re  g(z)  -  0  for  |z|  >  1,  we  find 


I(w)  -  0  for  |u-w  I  >  w_  (7) 

'  1  o 1  D 

Furthermore,  from  Eq.  (6)  it  follows  that  the  line  must  be  situated  at  the  zero 

A 

of  the  real  part  of  the  denominator,  e.g.,  at  the  solution  of  u  *  alm(g(w-coo)  + 

A 

g*(wQ-w)).  We  know  that  g(z)  attains  its  extreme  at  z  —  ±  1  (for  a  Debye  model 
we  have  g(±l)  —  ±®>)  and  that  w  >  »  a,  which  gives  for  the  position  of  the 

line  w  _  Together  with  Eq.  (7),  this  explians  the  sharp  edge  of  the 

line  at  the  blue  side.  In  conclusion,  the  finite  memory- time  in  the  atom-crystal 
interaction  prohibits  the  factorization  of  the  density  operator  in  thermal 
equilibrium,  which  in  turn  gives  rise  to  a  new  spectral  line  at  l>q  +  provided 

that  the  system  has  a  permanent  dipole  moment.  Its  profile  is  given  by  Eq.  (6). 

As  a  last  remark  we  mention  that  the  detailed  shape  of  the  line  at 
depends  crucially  on  the  adopted  model  for  the  phonon  dispersion  relation.  The 
sharp  peak  at  u  -  is  a  consequence  of  the  cutoff  of  the  dispersion 

relation  at  In  a  more  realistic  model,  where  the  detailed  structure  of  the 

first  Brillouin  zone  for  a  specific  crystal  is  taken  into  consideration,  this 
peak  will  smoothen  out.  Expression  (6)  for  the  line  shape  remains  valid,  but  the 
function  g(z)  will  behave  more  realistic  .  The  infinities  in  z  *  ±  1  will  turn 
into  finite  van  Hove  singularities.^ 


This  research  was  supported  by  the  Office  of  Naval  Research  and  the  Air 
Force  Office  of  Scientific  Research  (AFSC) ,  United  States  Air  Force,  under 
Contract  No.  F49620-86-C-0009 . 


TO 


/  A  V  A  A  A..1.  A  A. 


s’S'X'V's 


REFERENCES 

1.  Z.  W.  Gortel,  H.  J.  Kreuzer,  P.  Piercy  and  R.  Teshima,  Phys .  Rev.  B  7J_,  5066 
(1983) . 

2.  A.  A.  Maradudin,  E.  W.  Montroll,  G.  H.  Weiss  and  I.  P.  Ipatova,  Theory, of 
Lattice  Dynamics  in  the  Harmonic  Approximation.  Suppl.  3  of  Solid  State 
Physics  (Academic,  New  York,  1971). 

3.  S.  Efrima,  C.  Jedrzejek,  K.  F.  Freed,  E.  Hood  and  H.  Metiu,  J.  Chem.  Phys. 
79,  2436  (1983). 

4.  H.  F.  Arnoldus,  S.  van  Smaalen  and  T.  F.  George,  Adv.  Chem.  Phys.  (1987),  in 
press . 

5.  W.  H.  Louisell,  Quantum  Statistical  Properties  of  Radiation  (Wiley,  New 
York,  1973),  Ch.  6. 

6.  H.  F.  Arnoldus  and  T.  F.  George,  J.  Math.  Phys.  28,  2731  (1987). 

7.  H.  F.  Arnoldus  and  T.  F.  George,  Phys.  Rev.  B  (1987),  submitted. 

8.  A.  Pe remans ,  J.  Darville,  J.  M.  Gilles  and  T.  F.  George,  Phys.  Rev.  B  35, 
2690  (1987). 

9.  H.  F.  Arnoldus  and  T.  F.  George,  Phys.  Rev.  B  .36,  1  (1987). 

10.  A  low- intensity  laser  can  only  cause  transitions  between  two  states  |n>  and 

|m>  for  which  <n|H^|m>  *  0,  according  to  the  Golden  Rule.  The  Hamiltonian 
H^  from  Eq.  (1)  connects  only  states  which  differ  by  one  phonon.  See  for 
instance:  S.  W.  Lovesey,  Condensed  Matter  Physics.  Dynamic  Correlations 

(Benjamin,  London,  1980)  pg.  3. 

11.  N.  W.  Ashcroft  and  N.  D.  Mermin,  Solid  State  Physics  (Holt,  Rinehart  and 
Winston,  New  York,  1976)  pg.  465. 


Fig.  1.  Absorption  profile  I(w)  from  Eq.  (5)  as  a  function  of  w  -  and  for 

2 

u>o  -  3u)p  (dotted  line),  a  -  0.4  x  w^,  7  -  10  (low  temperature),  and 
0.5  x  The  left  peak  at  w  <  is  the  far  red  wing  of  the  regression  part  of 

the  line,  which  is  brought  about  by  transitions  with  w  -  wphonon  (diagram  (a)  of 
Fig.  2),  and  governed  by  the  transition  dipole  moment  The  Pea^  on  c^e 

right-hand  side  is  the  memory- time- induced  extra  resonance,  with  |w-a>o|  - 
^phonon'  Processes  represented  by  the  diagrams  (b)  and  (c)  of  Pig-  2  give 

rise  to  this  line,  and  its  line  shape  is  given  by  Eq.  (6).  In  a  zero  memory-time 
approximation  this  line  vanishes. 


Fig.  2.  Energy-conserving  diagrams  which  contribute  to  the  absorption  spectrum 
I(u>).  Double  arrows  indicate  a  photon  and  single  arrows  a  phonor..  Diagram  (a) 
and  its  reverse  are  responsible  for  the  profile  at  u>  <  whereas  diagrams  (b) 
and  (c)  represent  processes  which  yield  the  line  at  u>  >  Since  level  1 2>  is 
not  populated,  the  reverse  processes  (stimulated  emission)  of  diagrams  (b)  and 
(c)  do  not  occur. 
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